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Abstract 

Volatility modelling has become a significant area of research within Financial 
Mathematics. Wiener process driven stochastic volatility models have become popular 
due their consistency with theoretical arguments and empirical observations. However 
such models lack the ability to take into account long term and fundamental economic 
factors e.g. credit crunch. 

Regime switching models with mean reverting stochastic volatility are a new class 
of stochastic volatility models that capture both short and long term characteristics. 
We propose a new general method of pricing options for these new class of stochastic 
volatility models using Fouque's perturbation based option pricing method. 

Using empirical data, we compare our option pricing method to Black-Scholes and 
Fouque's standard option pricing method and show that our pricing method provides 
lower relative error compared to the other two methods. 
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1. Introduction and Outline 



Volatility modelling has become a significant area of research within Financial 
Mathematics: it helps us understand price dynamics since it is one of the key variables 
in a stochastic differential equation. Volatility is the only variable in the Black-Scholes 
option pricing equation that is unobservable, hence crucial to option pricing. Finally, 
volatility has a wide r ange of industria l applic ations from pricing exotic derivatives to 
asset pricing models RebO^ . Shiller Sch81 | also claims volatility can be used as a 



measure of market efficiency. 
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With volatility modelling, stochastic volatility has achieved growing importance as 
they capture a ric her set o f empirical and theoretical characteristics compared to other 
volatility models 



MR05bj : 



firstly, stochastic volatility models generate return distributions similar to what 
is empirically observed. For example, the return distribution has a fatter left tail 
and peakedness compared to normal distributions, with tail asymmetry controlled 



by p |Dur07 |. 



Secondly, Renault and Touzi [RT96l | proved volatility that is stochastic and p=0 
always produces implied volatilities that smile (note that volatility smiles do not 
necessarily imply volatility is stochastic). 

Thirdly, historic volatility shows significantly higher variability than would be ex- 
pected from local or time dependent volatility, which could be better explained by 
a stochastic process. A particular case in poi nt is th e dramatic change in volatil- 
ity during the 1987 October crash (Schwert Sch9cl | gives an empirical study on 



this). Finally, stochastic volatility accounts for the volatility's empirical depen- 
dence on the time scale measured, which should not occur under local or time 
dependent volatility. 

Wiener process driven stochastic volatility models capture price and volatility dy- 
namics more successfully compared to local and time dependent volatility models. 
However, for longer term dynamics and fundamental economic changes (e.g. "credit 
crunch"), no mechanism existed to address the change in volatility dynamics and it has 
been empirically shown that volatility is related to long term and fundamental condi- 
tions. To model volatility (and stock prices) in continuous time but also capture long 
term and fundamental factors requires regime switching models with mean reverting 
stochastic volatility. This is a currently growing area of stochastic volatility models. 

At present, there is no method for pricing options for volatility go verned by regime 



switching with mean reverting stochastic volatility. Using Fouque's |FPS00al | pertur- 
bation approach to stochastic volatility option pricing, we show how it is possible to 
obtain option prices. Additionally, the option pricing method is possible for any mean 
reverting process. We demonstrate our method with empirical results from S&P 500 
index options. 

The outline of the paper is as follows. Firstly we give a literature review of volatility 
models and the option pricing methods. We then explain Fouque's option pricing 
method under standard mean reverting stochastic volatility and then show how it 
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can be applied to regime switching with meaning reverting stochastic volatihty. We 
finally provide numerical experiments, demonstrating that our method of option pricing 
provides more accurate prices compared to standard Fouque option pricing (without 
regimes) . 



2. Literature Review on Volatility Models 

Bachelier BacOCll | proposed the first model for stock prices. Bachelier reasoned that 



investing was theoretically a "fair game" in the sense that statistically one could neither 
profit nor lose from it. Hence Bachelier inclu ded th e Wiener process to incorporate 
the random nature of stock prices. O sborne Osb59[ | conducted empirical work sup- 
porting Bachelier's model. Samuelson Sam65l | continued the constant volatility model 
under the Geometric Brownian Motion stock price model on the basis of economic 
justifications. 

Over time, empirical data and theoretical arguments found constant volatility to be 
inconsistent with the observed market behaviour (such as the leverage effect). A plot 
of the empirical daily volatility of the S&P 500 index clearly shows volatility is far from 
constant. This consequently led to the development of dynamic volatility modelling, 
which we will now discuss. 

2.1. Time Dependent Volatility Model 

It was empirically observed that implied volatility varied with an option's expiration 
date. Consequently, a straight forward extensi on prop osed to the constant volatility 
model was time dependent volatility modelling Wjl98|: 



dX/X = fxdt + a{t)dW. 



(1) 



Merton Mer73[ | was the first to propose a formula for pricing options under time 
dependent volatility. The option price associated with X is still calculated by the 
standard Black-Scholes formula except we set a = ac where: 



1 



T-t 



T 



a'^{T)dT, 



(2) 



i.e. di and c?2 in the Black-Scholes equation become: 



di 



do 



^o^?(f)+MT-t) + ^/fa2(r)rfr 
l09(^§^+f^{T-t)-^J^a\T)dT 



f^a'^{T)dT 



(3) 



(4) 



The equation ([2]) converts cr(t) to its constant volatility equivalent Oc over time period 
t to T. The distribution of X(t) is given by: 

logiX^lX,) ~ AT (^(/i - \al){T - t), al{T - t)^ . (5) 

Note that the constant volatility Cc changes in value as t and T change. This property 
enables time dependent volatility to account for empirically observed implied volatilities 
increasing with time (for a given strike). 

2.2. Local Volatility Model: Definition and Characteristics 

In explaining the empirical characteristics of volatility, a time dependent volatility 
model was found to be insufficient. For instance, time dependent volatility did not 
explain the volatility smile, nor the leverage effect, since volatility cannot vary with 
price. Therefore volatility as a function of price (and optionally time) was proposed, 
that is local volatility is a = f{X,t): 

dX/X = ndt + a(X, t)dW. (6) 

The term "local" arises from knowing volatility with certainty "locally" or when X is 
known -for a stochastic volatility model we never know the volatility with certainty. 

The advantages of local volatility models are that firstly, no additional (or un- 
tradable) source of randomness is introduced into the model. Hence the models are 
complete, unlike stochastic volatility models. It is theoretically possible to perfectly 
hedge contingent claims. Secondly, local volatility models can be also calibrated to per- 
fectly fit empirically observed implied volatility surfaces, enabling consistent pricing of 



the derivatives (an example is given in Dup97l |). Thirdly, the local volatility model is 



able to account for a greater degree of empirical observations and theoretical arguments 
on volatility than time dependent volatility (for instance the leverage effect). 

Some common local volatility m odels a re: the Constant Elasticity of Variance Model 
(CEV) pr oposed by Cox and Ross |CR76| . Mixture Distribution Model s by Brigo and 



Mercurio 



Kani |DK94 |. 



BMOOl ] and the Implied Local Volatility by Dupire Dup94| . Derman and 



2.3. Stochastic Volatility Model: Definition and Characteristics 

Although local volatility models were an improvement on time dependent volatil- 
ity, they possessed certain undesirable properties. For example, volatility is perfectly 
correlated (positively or negatively) with stock price yet empirical observations suggest 
no perfect correlation exists. Stock prices empirically exhibit volatility clustering but 
under local volatility this does not necessarily occur. Consequently after local volatility 
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development, models were proposed that allowed volatility to be governed by its own 
stochastic process. We now define stochastic volatility. 

Definition 1. Assume X follows the stochastic differential equation 



dX/X = fxdt + a{uj)dWi. 



(7) 



Volatility is stochastic if is governed by a stochastic process that is driven by an- 
other (but possibly correlated) random process, typically another Wiener process dW2- 
The probability space (fi, JF, P) is Q = C([0, oo) : M^), with filtration {J^t}t>o represent- 
ing information on two Wiener processes {Wi,W2)- 

The process governing cr{uj) must always be positive for all values since volatility can 
only be positive. The Wiener processes have instantaneous correlation p G [—1, 1] de- 
fined by: 

corr{dWi{t),dW2{t)) = pdt. 

Empirically p tends to be negative in equity markets due to the leverage effect but 
close to in the curr ency m arkets. Signi ficant s tochastic volatility r nodels i nclude: 
Johnson and Shaiino_ JS87| , Scott Mode]_|Sco87|, Hull- White Model |hW87| . Stein 
and Stein Model |SS9l|, Heston model |Hes93l | 

The key difference between local and stochastic volatility is that local volatility is 
not driven by a random process of its own; there exists only one source of randomness 
{dWi). In stochastic volatility models, volatility has its own source of randomness 
{dW2) making volatility intrinsically stochastic. We can therefore never definitely de- 
termine the volatility's value, unlike in local volatility. 

The disadvantages of stochastic volatility are firstly that these models introduce 
a non-tradable source of randomness, hence the market is no longer complete and 
we can no longer uniquely price options or perfectly hedge. Therefore the practical 
applications of stochastic volatility are limited. Secondly stochastic volatility models 
tend to be analytically less tractable. In fact, it is common for stochastic volatility 
models to have no closed form solutions for option prices. Consequently option prices 
can only be calculated by simulation. 

Stochastic volatility models fundamentally differ according to their driving mech- 
anisms for their volatility process. Different driving mechanisms maybe favoured due 
to their tractability, theoretical or empirical appeal and we can categorise stochastic 
volatility models according to them. Many stochastic volatility models favour a mean 
reverting driving process. A mean reverting stochastic volatility process is of the form 
FPSOoJ: 



= f{Y), 
dY = a{m-Y)dt + (3dW2, 



(9) 
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where: 



• P > and (3 is a constant; 

• m is the long run mean value of a; 

• a is the rate of mean reversion. 

Mean reversion is the tendency for a process to revert around its long run mean value. 
We can economically account for the existence of mean reversion through the cob- 
web th eorem, which claims prices mean revert due to lags in supply and demand 
LH92l |. The inclusion of mean reversion (a) within volatility is important, in particu- 
lar, it controls the degree of volatility clustering (burstiness) if all other parameters are 
unchanged. Volatility clusteri ng is a n important empirical characteristic of many eco- 



nomic or financial time series |Eng82| | , which neither local nor time dependent volatility 
models necessarily capture. Additionally, a high ^ can be thought of as the time re- 
quired to decorrelate or "forget" its previous value. 

The equation © is an Ornstein-Uhlenbeck process in Y. It is known the solution 
to equation ([9]) is: 

Y{t) =m+ {Y{0) -m)e-"' + (3 [ e'^^dWa, (10) 

Jo 

where Y(t) has the distribution 

Y{t) r^^^(m + (F(0) - m)e-"*, ^(1 - e'^"*)^ . (11) 

Note that alternative processes to equation could have been proposed to define 
volatility as a mean reverting stochastic volatility model, for example the Feller or 



Cox-IngersoU-Ross (CIR) process jFPSOObj : 



dV = a{m-Y)dt + pVYdW2. (12) 

However with a = f{Y) in equation we can represent a broad range of mean 
reverting stochastic volatility models in terms of a function of Y. 

2.4- Regime Switching Volatility 
2.4.1. General Regime Switching 



A class of models that address fundamental and long term volatilit y mode 



regime switching model (or hidden Markov model) e.g. as discussed in TimOOj, EvdH97 |. 



l ing is the 



In fact, Schwert suggests in |Sch89l | that volatility changes during the Great Depression 



can be accounted for by a regime change such as in Hamilton's regime switching model 
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Ham89| . For regime switching models, generally the return distribution rather than 
the continuous tir ne proc ess is specified. A typical example of a regime switching model 
is Hardy's model |HarOl[ |: 



logiiXit + l)/X(t))K) ~ Af{u,, cpi),te {1, .., i?}, 



(13) 



where 



• and Ui are constant for the duration of the regime; 

• i denotes the current regime (also called the Markov state or hidden Markov 
state); 

• R denotes the total number of regimes; 

• a transition matrix A is specified. 

For Hardy's model the regime changes discretely in monthly time steps but stochasti- 
cally, according to a Markov process. 

Due to the ability of regime switching models to capture long term and fundamental 
changes, regime switching models are primarily focussed on modelling the long term 
behaviour, rather than the continuous time dynamics. Therefore regime switching 
models switch regimes over time periods of months, rather than switching in continuous 
time. Examples of regime swit ching mo dels that model dynamics over shorter time 
periods are Valls-Pereira et al. VPHS04I|. w ho propose a regime switching GARCH 
process, while Hamilton and Susmel HS94l | give a regime switching ARCH process. 
Note that economic variables other than stock returns, such as infiation, can also be 
modelled using regime switching models. 

The theory of Markov models (MM) and Hidden Markov models (HMM) are meth- 
ods of mathematically modelling time varying dynamics of certain statistical processes, 
requiring a weak set of assumptions yet allow us to deduce a significant number of prop- 
erties. MM and HMM model a stochastic process (or any system) as a set of states 
with each state possessing a set of signals o r observ ations. The model s have been used 
in diverse applications such as econ omics |SSS02|| . queuing theor y |SF06l| . engineer- 



ing 



TG01| and biological modelling [MGPGof Following Taylor [TK84l | we define a 



Markov model: 

Definition 2. A Markov model is a stochastic process X{t) with a countable set of 
states and possesses the Markov property: 



Pilt+i = J I qi, q2, ..,qt = i) = p{qt+i = j \ qt 



(14) 



where 
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• qt is the Markov state (or regime) at time t of X{t); 

• i and j are specific Markov states. 

As time passes the process may remain or change to another state (known as state 
transition) . The state transition probabihty matrix (also known as the transition kernel 
or stochastic matrix) A, with elements aij, tells us the probability of the process 
changing to state j given that we are now in state i, that is aij = p{qt+i = j \ Qt = i)- 
Note that subject to the standard probability constraints: 

0<a,,- < l,Vi,i, (15) 

oo 

= l,Vi. (16) 

i=i 

We assume that all probabilities are stationary in time. From the definition of a MM 
the following proposition follows: 

Proposition 1. A Markov model is completely defined once the following parameters 
are known: 

• R, the total number of regimes or (hidden) states; 

• state transition probability matrix A of size RxR. Each element is aij = p{qt+i = 
j\qt = i), where i refers to the matrix row number and j to the column number of 
A; 

• initial (t=l) state probabilities tTj = p{qi = i),\/i. 

A hidden Markov model is simply a Markov model where we assume (as a modeller) 
we do not observe the Markov states. Instead of observing the Markov states (as 
in standard Markov models) we detect observations or time series data where each 
observation is assumed to be a function of the hidden Markov state, thus enabling 
statistical inferences about the HMM. Note that in a HMM it is the states which must 
be governed by a Markov process, not the observations and throughout wc will assume 
one observation occurs after one state transition. 

Proposition 2. A hidden Markov model is fully defined when the parameter set {A, B, tt} 

are known: 

• R, the total number of (hidden) states or regimes; 

• A, the (hidden) state transition matrix of size RxR. Each element is a^j — 
P((lt+i ^j\qt^ i); 

• initial (t=l) state probabilities TTj = p{qi = i),\/i; 
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• B, the observation matrix, where each entry is bj{Ot) = p{Ot\j) for observation 
Ot- For bj{Ot) is typically defined to follow some continuous distribution e.g. 
hj{Ot) ~ N'{uj,Lpj). 

Regi me swit ching; has been developed by various researchers. For example, Kim 
and Yoo KY95l | de velop a multivariate regime switching model for coincident economic 
indicators. Honda Hon03l | determines the optimal portfolio choice in terms of utility 
for assets following G BM bu t with continuous time regime switching mean returns. 
Alexander and Kaeck AK08I apply regime switching to credit default swap spreads, 



Durland and McCurdy DM94l | propose a model with a transition matrix that specifies 
state durations. 

2.4.2. Regime Switching Mean Reverting Stochastic Volatility Models 

A common shortfall in all the volatility models reviewed so far has been that all 
these models are short or long term volatility models. The models implicitly ignore 
any long term or broader economic factors influencing the volatility model, which is 
empirically unrealistic and theoretically inconsistent. 

Combining Wiener process driven stochastic volatility models with regime switching 
would give the benefits of capturing the short term price dynamics while the regime 
switching would capture volatility changes due to fundamental and longer term effects 
in a tractable method. In fact Alexander combines two models to model short and 



long term dynamics, giving the binomial mixture distribution model Ale04l | (mixture 
distribution model and a binomial tree). 

A model that addresses both the long and short term dynamics are regime switching 
mean reverting volatility models (RSMR). These RSMR models are of the form where 
the intra-regime (within regime) dynamics follow an exponential Ornstein-Uhlenbeck 
process (expOU): 



dX/X 
da'Ja'i 



Hidt + (j[{uj)dW] 



1; 



dt + (3[dW2, 



(17) 
(18) 



I e {1,2}. (19) 

Alternatively, it can be more convenient to express RSexpOU in terms of OU: 

dX/X = fxdt + ai{uj)dWi, (20) 
ai = exp{Yi{t)), (21) 
dYi{t) = ai{mi - Yi)dt + f3idW2. (22) 

We denote i as the current regime and for convenience set R=2. The regime switch- 
ing process is a discrete time Markov process (e.g. changing in time intervals of one 
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year). We note that regime switching models that change state continuously do not 
capture long term and sustained changes. In fact Hamilton's regime s witchin g models 
have discrete time switching periods of three months (see for instance EH90|). 

The literature on regime switching models with stochastic differential equations 
in their continuous time price dynamics is currently growing but is not signi ficant , 
especially f or mean reverting stochastic volatility. Examples include Kalimipalli KS04l | 
and Smith Smi02|, who describe a st ochastic volatility model with regime switching 



to model interest rates, Elliot et al. ECS05| prop ose a regime switching Geometric 
Brownian Motion model (GBM) and Mari MarOGl ] models electricity spot prices with 
each regime possessing a different stochastic differential equation structure, such as 
mean reversion in the drift. 



3. Perturbation Based Option Pricing 

We price options for RSMR model using Fouque's option pricing method, which is 
based on perturbation theory. Fouque's perturbation method is normally applied to 
mean reverting stochastic volatility models; we briefly introduce Fouque's perturbation 
theory by flrst introducing perturbation analysis and then Fouque's method. Finally 



we show we can apply Fouque's method to RSMR mode 



s. For more information on 



perturbation theory the reader is referred to Hinch Hin91|, Simmonds SM98l |. Bender 
and Orszag |b099 |. 



3.1. Introduction to Perturbation Analysis 

For many mathematical equations we cannot flnd precise solutions, however, we 
would like to obtain approximations in a mathematically consistent and rigorous way. 
The perturbation method is applied when we know the solution to an equation but 
would like to know the solutions when variables are minutely changed (or perturbed) 
by a small amount e, where e ^ but is flnite. Perturbation problems can be solved by 
an iterative method or by a series expansion method; Fouque adopts the latter method. 

Perturbation analysis can be explained by example; consider the following example 
equation we wish to solve: 



X 



X. 



(23) 



This equation can be easily solved using knowledge of quadratics. Now consider when 
we perturb this problem by e: 



X 



ex, where e ^ 0. 



(24) 
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Solving for x now becomes difficult for various small (but finite) values of e. Note that 
in this example it is still possible to find the 2 analytic solutions to equation using 
our knowledge of quadratics: 



Xa = e/2 + v/(l + e2)/4, (25) 
Xh = e/2 - v/(l + e2)/4. (26) 

However, instead of using knowledge of quadratics we apply the perturbation method: 
we determine the solutions Xa, x^ by approximating x by a series expansion. One could 
use any series expansion (e.g. Taylor, Maclaurin) but we use a power series expansion: 



oo 

X = ^x„e" (27) 

n=0 

xo + exi + 6^X2 + e^X3 + .... (28) 



where Xi,X2,X3 are constants to be determined. Note that xq is called the leading 
term, Xi the first correction, X2 the second correction and so on. The unperturbed or 
known solution corresponds to e = and therefore equals xq. Normally obtaining a 
solution to the unperturbed problem (that is xq) is easy to obtain, otherwise the set of 
equations (as shown later) that need to be solved to obtain approximations to Xa,Xfe 
will be difficult to solve. 

With the series approximation of x (equation (l28l) ). we insert this approximation 
into the original equation (IMll : 

x^ — 1 = ex 

(xo + exi + e^X2 + e^aJa)^ - 1 = e(xo + exi + t^X2 + 6^X3) 
(x^ + 2exoXi + e2(x? + 2x0X2) + 0(e3) + 0(e^))-l = e(xo + exi + ^(e^)). 

Now it is standard practice to ignore terms of order (denoted 0{e')) and higher 
since e is small and so terms of 0{e^) or higher are considered insignificant. Therefore 
the equation simplifies to 

Xq + 2exoXi + e2(x^ + 2x0X2) — 1 = exo + e^Xi (29) 
(xg — 1) + e(2xoXi — Xo) + e2(x^ + 2x0X2 — Xi) = 0. (30) 

Now equation ( l30l) is significant because the bracketed terms represent coefficients of 
e for each power and each bracketed term separately must equal for the equation to 
balance. So now we can equate terms of the same order from equation flHUj) : 

0(e°) : x^ - 1 = 0, (31) 
O(e^) : 2xoXi-Xo = 0, (32) 
0{t^) : xl + 2x0X2 - xi = 0. (33) 
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Solving the equations (!3Ti) - (l33l) and remembering there are 2 solutions gives: 



Xa : {a;o = l,Xi = l/2,X2 = l/8}, (34) 

Xb ■ {xo = = 1/2,X2 = -1/8}. (35) 

Using these values, we can now calculate approximate but extremely accurate solutions 

to Xai Xb- 

X = Xq + exi + e^X2 + ■■■■ (36) 

Xa = l + e(l/2) + 6^(1/8) (37) 

Xb = -l + e(l/2) + e2(-l/8). (38) 

To summarise, the three steps of perturbative analysis (by Bender and Orszag 



B099| ) are: firstly to convert the original (unperturbed) problem into a perturbed 
problem by introducing e. Secondly, we assume the expression for the solution is in 
the form of a perturbation series and determine the coefficients of that series. Thirdly, 
we determine answers to the solutions by summing the perturbation series for the 
appropriate values of e. 

Perturbation problems can be divided into two categories: regular and singular 
perturbation problems. Singular pertubation problems have significantly different so- 
lutions as e — > when compared to the unperturbed problem (e = 0). If the solution 
changes insignificantly with e then we have a regular perturbation; the perturbed so- 
lution smoothly approaches the unperturbed solution as e — 0. Fouque's perturbation 
based option pricing method is a singular perturbation problem. 

To give an example of perturbation problems, a spinning drawing pin rotates around 
its vertical axis and over a short time period its horizontal distance from its vertical 
axis is unlikely to be significantly different from its initial position (at t=0). However 
over long time periods, the drawing pin will eventually come to a physical rest and its 
horizontal position will be appreciably different than from at t=0. Therefore over a 
small time period the position can be modelled as a regular perturbation problem but 
over long time periods it is a singular perturbation problem. 

3.2. Perturbation Formulation of Option Pricing 

Option pricing where the stock p rice proce ss is governed by stochastic volatility can 



be analytically challenging. Fouque FPSOOal ] finds a method of pricing options whose 



underlying is governed by a mean reverting stochastic volatility process by solving 
the option pricing partial differential equation as a singular perturbation problem. 
A significant benefit of Fouque's method is that the option pricing method is model 
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independent; it can be applied to any fast mean reverting stochastic volatility model 
and without explicit knowledge of the volatility dynamics. Additionally, the option 
pricing calibration is far more parsimonious compared to other stochastic volatility 
option pricing methods. 

To find the option pricing partial differential equation when the underlying asset is 
driven by mean reverting stochastic volatility, we can use a replicating portfolio argu- 
ment. Alternatively, we can determine the option pricing partial differential equation 
associated with an option by risk neutral valuation and apply the Feynman-Kac equa- 
tion, which is the approach we take here. We abbreviate call option C(X,K,s,T,r,f(Y)) 
to C(s,X,Y), which is an option at time s with volatility as a function of some mean re- 
verting stochastic volatility process (f(Y)). By risk neutral valuation C(s,X,Y) is given 
by: 

C{s, X, Y) = e-'^^-'^E'^[X{T) - K| J^,]+. (39) 

We recall that in a no arbitrage market under stochastic volatility there is more than 
one risk neutral measure Q (no unique Q measure exists). The possible Q measures can 
be parameterised by the market price of volatility risk 7 = 7(t, X, f{Y)). Therefore Q 
measures are a function of 7, consequently the price of option C(t,X,Y) is a function 
of 7. 

To obtain the (risk neutral) option pricing partial differential equation by the 
Feynman-Kac equation, we require the risk neutral process. We begin with a stock 
price with a mean reverting stochastic volatility process, that is: 

dX/X = fxdt + adWi, (40) 
^ = f{Y), (41) 
dY = a{m-Y)dt + (3dW2. (42) 

The risk neutral process for dX/X is: 

dX/X = rdt + adW^. (43) 



We can re-express dW2 as |FPS00a|: 



dW2 = pdWi + Vl - p'^dW*, (44) 

where W* is a standard Wiener process independent of Wi. To change dY to its risk 
neutral process we require: 
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1. 



2. 



i.e. dWi 



i.e. dW* 



= dWi + ^^dt, 
- dW^-^dt. 

dW* + -f{t,X,Y)dt, 
--f{t,X,Y)dt. 



To obtain the risk neutral process of dY we substitute in equations 
rearrange dY; starting with substituting equation (jSj) for dW2 gives 



dY 



a{m - Y)dt + pdW2, 

Y)dt + (3 \pdWi + - p'^dW* 



aim 



Substituting equation ( l46l) for dWi and equation ( HSl) for dW* gives: 
dY = a{m - Y)dt 



+ /5 

Rearranging now gives: 
dY 

+ /5 



-Y)-(3 [py^ + V^^lit, X, Y) 



(45) 
(46) 



(47) 
(48) 

HI) and 



Since the risk neutral equation for dX is a two dimensional diffusion equation in 
iyV'^^W'^*), the associated option pricing partial differential equation for C(s,X,Y) is 
obtained by the application of the multidimensional version of the Feynman-Kac equa- 
tion with boundary condition C(T,x,y). 



The Feyman-Kac equations give the option pricing partial differential equation: 




9C ^1.2^ ^ 2d'C 



+ pPxfiy)-^ + + r (x— - C 

dxdy 2 dy^ \ dx ^ 



+ [a{m-y) - pA{s,x,y)] 



9C 
dy' 



where A(s, x, y) 



P^-j^ + V^'l{s,x,y). 
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The term A{t,x,y) can be interpretted as a combination of market price of risk and 
market price of volatility risk. Fouque assumes that A(s, x, y) = A{y), a function of y 
only: 

My) = P^-^^ + Vi^Hy)- (49) 
f{y) 



the distribution becomes Af{m, |-) (see equation (fTTl)). If we denote the variance of 



Now from our knowledge of OU processes (which is also Y(t)) we know as t — oo 

^^ 

this long term distribution as f^, that is t;^ = then (3 = vy/2a and we can re-express 
the option pricing partial differential equation as: 

= -^ + -r{y)x'j-^ + pv^2axf{y)—— + v^a—^ + r[x—-C 
OS 2 ax^ oxoy oy^ \ ox 

dC 

+ [a{m-y) -vV2aA{y)]-^. 

Now if we re-express a as a = -, where e is small but finite (since a oo) 
our equation becomes: 

dC , 1 ,2, N , vV2 , d^C , v^&'C , f dC 



ds ^ 2^ ^^^^ dx"^ ~^ ^ y/e ^^^^^ dxdy ~^ e dy"^ ^ ^ \ dx ^ 



dC_ 
dy ' 



We can also re-express the equation using the following operator notation for con- 
venience: 

Co = -"£-. + (^-y)^^ (50) 

A = V2pvxm^-V2vA{y)-^, (51) 

Therefore the option pricing partial differential equation is: 

-Co + ^Ci + C2]c = 0. (53) 
/ 

Equation fl53l) is a (singular) perturbation problem; we can therefore solve the option 
pricing partial differential equation (under mean reverting stochastic volatility) by 
perturbation methods. We note that: 

• we sometimes express C2 = CBsifiy)), where Cbs{-) is the Black-Scholes partial 
differential operator with volatility (■). 
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• the Co operator gives a Poisson equation for Y (an OU process). Note that the 
popular definition of a Poisson equation is a form of Laplace equation (second 
order partial differential equation), however, for stochastic processes it is Cq. 

3.3. Perturbation Solution to Option Pricing 

To determine the solution to C(s,x,y) by a perturbation approach we expand C in 
powers of y/e 

C = Co + v^Ci + tC2 + eiCg + (54) 

and insert equation into equation (!53l) : 

= Q^o + ^^Ci + a) (Co + v^Ci + eC2 + eiC3 + ) 

= —CqCq 

e 

+ ^(^oCi + ACo) 

+ (£0^2 + ACi + ACo) 
+ ^t{C,C, + C,C2 + C2C,) 
+ 0{e). 

As with standard perturbation problems, we can equate terms of the same order to 
assist us in solving the equation. Fouque solves C upto C = Co + \/eCi (although we 
could solve further), therefore our aim is to find Co and Ci. 

Zero Order Term (Cq) 

We begin by firstly determining Co; by equating terms of order 1/e we have: 

>CoCo = 0. (55) 

For equation ( l55ll to be correct and since Cq acts only on y, Co must be a function of 
s,x only: Cq = Co(s,a;). 

Now equating terms of order 1 / ^/e we have: 

CqCi + CiCq = 0. (56) 

For equation (|56ll to be correct, if Ci takes derivatives on y only and since we have 
already deduced Co = Co(s,a;) therefore we have CiCq = 0. This in turn implies from 
equation fl3Ul) that CqCi = 0, which in turn implies Ci = Ci(s, x) for equation fl3U|) 
to be correct. If desired one could continue to apply this iterative method for higher 
orders 1, e and so on to obtain further equations. 
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Now the order 1 terms give: 

£oC2 + CiCi + C2C0 = 0. (57) 

Since we have aheady reasoned CiCi=0 we can write 

C0C2 + C2C0 = 0. (58) 

Regarding the variable x as fixed C2C0 is a function of y since C2 contains f(y); if we 
focus on the y dependency only we can rewrite equation fl3H|) as: 



Coxiy) + g{y) = 0, (59) 

where g{y) = C2C0. Fouque states this is a Poisson equation in x{y) and so can assert 
for a solution to exist g(y) must fulfil the "centering condition": 

POD 

<g{y)>= / giyMy)dy = o. (60) 



The terms in equation (1601) relate to the theoretical area of ergodicity and will be 
covered in more detail in section 13. 4t we briefiy define them for now. We define < 
g{X) > as: 



<g{X) >= Mm I f g{X)ds. 



(61) 



For any process X that is ergodic or a function g() of an ergodic process g(X), there 
exists an invariant distribution or long term distribution exhibited by g(X) as time 
tends to a large limit. We denote $(■) as the probability measure of the invariant 
distribution. Therefore 

/oo 
g{X)^{X)dX. (62) 
-00 

Equation ( I6OI) therefore means < g{y) > equals its average under its invariant distri- 
bution; note that both < giij) > and E^[g{X)\ are constants. We also denote 

{af = <f{Y)> (63) 

when f{Y) = a, that is when we operate <> upon the squared mean reverting stochas- 
tic volatility process. The vari able a is a constant and is referred to as the "effective 



volatility" by Fouque |FPS00a|. 



From equation ( l59l) Cox = ~g{y) therefore substituting this into equation ( I6OI) 
gives: 

<C2Co> = 0. (64) 
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Since we have already established Cq does not depend on y we have 

<£2>Co = 0. (65) 

Now =< P{y) > and so by definition of Cbs we have 

<C2> = Cbs{o), (66) 
so that Cbs{^)Cq = 0. (67) 

Therefore we finally deduce the value of Co'- it is the option price obtained by the 
Black-Scholes equation with constant volatility a = a. 

First Order Term (Ci) 

Since the centering condition < C2C0 >=0 is satisfied we can re-express C2C0. Recall- 
ing that we can denote C2C0 = CBs{fiy))Co we have: 

C2C0 = C2C0- < C2C0 >, since < C2C0 >= 0, (68) 
= CBsifiy))Co-CBsia)Co (69) 

From equation 0581) we have: 

C0C2 + C2C0 = (71) 
C0C2 = -C2C0. (72) 



Now substituting C2C0 expression from equation (ITOl) we have an expression for C2: 

= -r^Co (73) 
lim-^V-'-^) (74) 

so that = ('!(/%) _jf=)x2^^') (75) 

= -i£„-'(f(,)-aV^. (76) 
The last line was possible since jCq takes derivatives in y only. If we write 

CocP{y) = f{y)-a\ (77) 

then we have a Poisson equation for which (f){y) is a solution to it. We can therefore 
express equation (17H|) as 

C2 = -l{<P{y) + k{s,x))x'^, (78) 



where k(s,x) is some constant with respect to y, that may depend on (s,x). 
Now recall that we have: 

£oC3 + + ACi = (79) 
sothat^oCs = -(£1^2 + £2^1). (80) 

Equation (IHOj) is another Poisson equation but in C3, consequently we can apply the 
centering condition: 

< £1^2 + £2^^! > = (81) 

< £2C'i > = — < £16*2 >, since <> is an integral. (82) 

Now 

1. Ci is not a function of y and it was aheady mentioned < £2 >= C-bs^) so we 
have 

< £2^1 > = < £2 > a (83) 

= Cbs{^)Ci. (84) 



2. substituting our expression for C2 from equation (ITHl) . we have: 

- < £iC2 >=\< > x'^. (85) 

Note that the term k(s,x) disappears since £1 takes derivatives on y terms only. 
Using equation (|82l) we can equate — < £iC2 > and < C2C1 >: 

Cbs{^)C^ = ^<C^<P{y)>x'^. (86) 
Now recalling the Ci operator definition (equation (CTj) we can determine < Ci(p{y) >: 

V2pvf{y)x-^ - V2vA{y)—\ 0(y) > (87) 
< V2pvf{y)(f)'{y)x-^ - V2vA{y)(j)'{y) > (88) 
V2pv < f{y)<P\y) >x-^-V2v< A{y)<P\y) > . (89) 
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We insert equation (18911 into equation (18611 giving: 
£ss(?)Ci 



1 o'^C 

-V2pv < f{y)<P\y) > + V2pv < f{y)<P'iy) 



> X' 



dx"^ 



V2 



V < A{y)(t)'{y) > X- 



2 " ~ ' " dx^ 

= -^pv<f{y)<p'{y)>x'^ 

+ (^V2pv < f{y)ct>'{y) > -^v < A{y)<P'{y) > j x'^. 



Recalling a = 1/e we can re-express Cbs{<^){V^Ci) in a more convenient form: 

Css{a){^eC,) = ^^-^^ + ^3-^^, (90) 

where 

V2 = -^{2p<fiy)<P'iy)>-<A<P'iy)>), (91) 
V2a 

Vs = -^<fiy)<P'iy)>, (92) 



2a 

with term inal condition Ci{T,x) = 0. It can be shown the solution to -^/eCi gives 
FPSOoJ: 



Therefore: 



C 



Co + V^Ci + eC2 + .... 
Co + v^Ci 



Co - (T - s) V2X 



(93) 

(94) 
(95) 
(96) 



FPSOOa 



As stated by Fouque 
fact is further explained below. 



dx"^ dx^ 
detailed expressions of V2, V3 are not important. This 



Option Pricing Using Observed Data 

The terms V2, V3 become explicit functions of the stochastic volatility if we fully define 
f in f(Y). However, Fouque proves that we can price C (specifically equation (lUUjl ) with- 
out specifying the stochastic volatility dynamics by using empirically observed data. 
To price C from equation ( l96ll we need to calculate Co, V2, V3, d'^Co/dx^ and d^Co/dx^. 
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The variable Cq is calculated using the standard Black-Scholes equation (where volatil- 
ity would be a), variables d'^Co/dx^ and d^Co/dx^ are known as the option's Gamma 
and Epsilon respectively and are calculated by 

^ = (97) 

dx^ xa^l-niT - s) 

and 

1 + , (98) 



dx^ x^a^lT^iT - s) V o^/T 

where d\ is taken from the Black-Scholes equation. 

The variables V2, V3 are found by fitting an affine function to a logarithmic plot of 
empirical option prices. To prove this, we first recall that implied volatility X is given 
by 

Cs5(x,s,T,r,J, i^) = C°''^ 

where Cbs is the European call option price under Black-Scholes option pricing. Now 
Qohs jg ]2iodelled by our perturbation expansion, that is: 

CBs{x,s,T,r,Z,K) = C"'' (99) 

= Co + v^Ci + .... (100) 

We can expand X as: 

I = Jo + v^Ji + ... (101) 

We can also use a Taylor series expansion to expand f unction Cbs- A function f(x) 



expanded around x=k by Taylor series expansion gives Kre93l | : 



Therefore a Taylor series expansion of Cbs{^i s-, T, r, I, K) around I = Iq gives 

CBs{x,s,T,r,T,K) = CBsix, s,T,r,Io, K) (103) 

da 

= CBs{x,s,T,r,Io,K) + y/eli — (105) 

+ .... (106) 

The last line was possible since y/e2i = I — Iq from equation (110 II) . 
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Equating equation (11051) with expansion C = Co + \ftC\ + .. we have: 



Cbs{x, s, T, r, To, K) + y^Ii 



da 



+ ... = Co + v^Ci + .... (107) 



By standard perturbation analysis we equate terms of the same order; matching terms 
of order 0{1) we have: 

CBs{x,s,T,r,Io,K) = Co^Io = a. 



Equating terms of order C(^/e) gives 



dCBs{x,s,T, r,a,K) 
da 



(108) 



If we now insert this expression for ^/eXl into our expansion for I we have: 



I = a + ^/^Ii + 0{e) 



a 



dCBs{x,s,T, r,,a,K) 
da 



+ 0{e). 



(109) 
(110) 



We can also make the substitution for ^^g^ in equation (1 11 01) using the option's Vega: 

-4 

OCbs 3;e 2 ^/T — s 



da 
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Now substituting equations (P7j) and (lUHj) into equation flU^ and after some rearranging, 
it can be shown we have: 



aL+b + 0{e) 



;il2) 



where 



log{K/x) 



a = — ^ ^ V3 = —aa 
a-^ 

h = a + ^(r+ '^a^\ - Yl^y^ = a^a - b) - a{r + fa^)). 



We can therefore obtain a, b from a simple linear plot of equation (11121) (given we know 
r and a). Plotting I on the y-axis and (log(K/x)/T) on the x-axis, the gradient gives 
a and the y-intercept is b (in fact we can think of b as the at the money X). 

We know that any function of Y f(Y) (that is any mean reverting stochastic volatility 
process) tends to a constant a"^ =< f^{Y) > and a can be measured from empirical 
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price data without specifying the volatih ty's dyna mics. Methods of measuring a are 



not required but the reader is referred to FPS00a| for more information 



We can now determine stochastic volatihty option prices under mean reversion 
without specifying the volatihty's dynamics. The variable a is obtained from price 
data, variables required to determine V2, V3 are observable (from the linear plot of 
equation flll2p ). Therefore Fouque's option pricing method is model independent; 
all that is required is that the stochastic volatility process is a function of Y (mean 
reverting). 

We note that the option pricing accuracy is of the order of 0{e): 

|C-(Co + v^Ci)| ^O(e). (113) 

Since ol = \ and we have assumed fast mean reversion we can assume the pricing error 
is negligible. For higher order approximations the option pricing method is no longer 
model independent and we would require specific knowledge of the stochastic volatility 
dynamics. 

3.J^. Regime Switching Perturbation Based Option Pricing 

Regime switching models generally specify the regime's distribution but not the 
intra-regime dynamics. Consequently if we can price intra-regime options without ex- 
plicit knowledge of the intra-regime dynamics (just assuming mean reverting stochastic 
volatility) it would be applicable to general RSMR models. 

Now Fouque's option pricing method is model independent as it does not require any 
specific definition of the mean reverting stochastic volatility process to price options. 
All the pricing variables are either exogenously determined variables (such as the strike 
and expiry) or are empirically observable (such as the stock price and interest rates), 
except 0", which is calculated from stock price data. Therefore if we can relate a to a 
typical regime switching model's specification we can obtain intra-regime option prices 
for general regime switching models, assuming a mean reverting stochastic volatility 
intra-regime process. We will now demonstrate this. 

If the intra-regime dynamics had simply been time dependent volatility (T{t), an 
analytic relation exists between a{t) and the return distribution. If we denote the 
constant volatility equivalent of cr{t) by cXc 



1 '•^ 



then the distribution would be given by: 

1 



log{X{T)/X{t)) -^f[[^^-^-al]{T-t), al{T -t)]. (115) 
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Our approach will be to demonstrate a similar relation for mean reverting stochastic 
volatility. 

In section 13.21 we defined a mean reverting stochastic volatility process where a is 
the rate of mean reversion and is empirically observed to be very high [a oo). A 
mean reverting stochastic volatility process a = f{Y) has the important property of 
ergodicity and it is this property which enables us to link a regime's distribution to its 
intra-regim e dynam ics. For a more detailed review of ergodicity the reader is referred 



to Kannan Kan79 | 



To explain ergodicity let g(X) be a process, where g is a function of an ergodic 
process X(t). The expectation for g(X) is 

/oo 
giX)piX)dX. (116) 
-oo 

Note that the expectation is generally a function of time. For an ergodic process X(t) 
or a function of an ergodic process g(X) there exists an invariant distribution, which 
is the long term or equilibrium distribution exhibited by g(X) as time tends to a large 
limit. If we denote $(■) as the probability measure of the invariant distribution for X(t) 
then the ensemble average (also known as the statistical average) is the expectation 
under $(■) probability measure: 

/oo 
g{XMX)dX. (117) 
-oo 

Now the time average of g{X) is defined as: 

W) = \ j\{r)dT. (118) 

Note that a time average is generally a random variable. We define the time average 
< g{X) > as: 

<g{X)>= hml fg{X)ds. (119) 

We say X(t) is an ergodic process if 

E''[giX)]=<giX)>, (120) 

that is X(t) is an ergodic process if g(X) has the property that the time average as 
t oo or < g{X) > equals its ensemble average under its probability measure $(X). 
Note that in equation fll20p both < g{X) > and E'^[g{X)] must always be a constant 
for the equation to hold. 
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Now for mean reverting stochastic volatility, a = f{Y) is ergodic therefore its time 
average approaches < fiY) > (and this is a constant): 



1 



hm - / f{Y)ds =< f{Y) > . (121) 

This relation is true regardless of the value of a. If t is large but finite and a oo 
then we have: 



' f{Y)ds^<f{Y)>. (122) 



Therefore the time average of a highly mean reverting stochastic volatility process 
approaches the constant < f(Y) >. 



The ergodic relation of equation fll2Up is also true for g'^{X) 

1 r 



FPSOOa 



< g\X) >= hm ^ [ g\X)ds. (123) 
Therefore the time average of cr^ = /"^{Y) approaches < g'^{Y) > (a constant): 

hmi f f{Y)ds = <f{Y)>. (124) 







We also defined in section 13.31 



t-+oo t 

If t is finite but a ^ oo then we have: 

1 



{af = <f{Y)>, (125) 
lim - / f{Y)ds. (126) 



'f\Y)ds ^ <f{Y)>, (127) 



(128) 

We also know that it is possible to convert time dependent volatility to its constant 
volatility equivalent Uc from section [2m that is equation (11141) . The constant volatility 
equivalent o"c for f(Y) is therefore approximated by: 



1 



= \lj I P{Y{s))ds, (129) 

a, ^ V< P{Y{s)) >, (130) 
^ a. (131) 

Therefore applying equation (11151) with cTc ~ a, we have a way relating any general 
mean reverting stochastic volatility process to the return distribution: a is the regime's 
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standard deviation. The distribution of each regime is approximated by equation f lllSp 
with (Tc ~ and as t ^ oo the approximation is precise. The variable a is specific to 
each regime since with each regime the volatihty parameter values may change (and 
therefore the return distribution); we therefore write a-i where i G {1, .., R}- 

Now as we are pricing intra-regime options we do not price options with expiries 
beyond the date of the next possible regime change. Consequently, for intra-regime 
options the regime does not change during the option's life and so we do not need to take 
into account risk arising from switching regimes. Therefore regimes do not introduce 
any additional incompleteness into intra-regime option pricing. Furthermore, a takes 
on only one value during the life of an option and so we can obtain intra-regime option 
prices using Fouque's option pricing method. 

Now as we have established Fouque's option pricing method is model independent, 
Wi (where i G {1, R}) is specified for general regime switching models and assuming 
we have RSMR, we can therefore price intra-regime options for any general RSMR 
model. Since we have 

C ~ Co + v^Ci + eC2 + ...., (132) 

our intra-regime option pricing equation for RSMR is therefore: 

Q = Co. + v^CH + eC2i + .... (133) 
~ Co^ + V~eCu,le{l,...,R}, (134) 

where a, is specified for each regime. 

If we wished to price options at any point in time without knowledge of the current 
state and beyond the next state transition then w e can apply Boyle's and Draviam's 



regime switching option pricing equation BDOTj], which is a transition probability 



weighted sum of option prices. However, long dated options (options with expiries 
longer than the duration of one regime (one year)) tend to be rarely traded, conse- 
quently their prices are significantly distorted by illiquidity effects (to be discussed in 
section |4?T!) . In fact many traded and liquid options expire at times considerably less 
than one year, hence the limit on option expiry is not restrictive. 



Option Pricing Advantages 

We now discuss the advantages of our option pricing method. Firstly, current regime 
switching option pricing methods tend to neglect intra-regime dynamics but also are 
model specific. Intra-regime option pricing using Fouque's perturbation approach is 
applicable to general regime switching models without requiring explicit knowledge of 
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intra-regime dynamics. All that is required is that we know (Xj and we assume mean 
reverting stochastic volatility intra-regime dynamics (since Fouque's option pricing 
method is model independent). 

A possible method of intra-regime option pricing would be to apply the appropriate 
option pricing formula to the intra-regime dynamics. However, this requires explicit 
knowledge of the intra-regime process and is not strictly regime based option pricing 
since such a method is not purely a function of the regime switching model. There 
is also the significant disadvantage that Wiener process driven stochastic volatility 
models with correlated Wiener processes have no analytic solutions, excluding the 
Heston model. Furthermore, the Heston model assumes a specific stochastic volatility 
process, risk neutral measure and there is no stated relation between the volatility of 
a regime switching model and the Heston model's volatility process. 

Secondly, our method assumes mean reverting stochastic volatility intra-regime dy- 
namics. The assumption of mean reverting stochastic volatility intra-regime dynamics 
is less restrictive and more realistic compared to local and time dependent volatility 
intra-regime dynamics. Hence such an intra-regime option pricing approach should pro- 
vide better option pricing compared to other intra-regime or general regime switching 
pricing methods. 

Thirdly, in contrast to other stochastic volatility option pricing methods or regime 
switching option pricing methods with no intra-regime dynamics, Fouque does not 
choose some particular function for specifying a unique risk neutral measure. Instead 
Fouque uses the "market's view" of the market price of volatility risk by plotting 
equation flll2p to extract a unique risk neutral measure from empirical option data. 
Consequently our intra-regime option pricing can be calibrated to current option data 
and is not restricted by any specific risk neutral measure. 

Fourthly, an important advantage of Fouque's option pricing method is that cali- 
bration is far more parsimonious compared to other stochastic volatility option pricing 
methods. All that is required is a and a linear fit of equation (11121) : the calibration is 
just as parsimonious under regime based Fouque option pricing where we use cTj instead. 
Typically Wiener process driven stochastic volatility option prices require calibrating 
many volatility parameters. 

Fifthly, using our knowledge of perturbation analysis we can use Fouque's method 
to obtain approximations to intra-regime options for RSMR. We can express option 
pricing for RSMR using equations (I133p and (11341) . From perturbation theory we know 
we can therefore approximate Cj by a high degree of accuracy by Coi since Coi is the 
largest term and the remaining terms are multiplied by e", where n > 0. Now we know 
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from equation (1671) that Co is simply the option price using the Black-Scholes equation 
therefore Cqj is simply the Black-Scholes option price for regime i with volatility a = o"j. 
This approximation has the advantage that we do not require option data for calibration 
or the market price of volatility risk. 

The approximation Ci ~ Cqi is advantageous since we can apply various regime 
switching option pricing formulae that ty pically a ssume volatility is constant in each 



regime, for instance Mamon and Rodrigo MR05al |. Additionally, C^i is calculated by 



the Black-Scholes equation, which is advantageous because much literature exists on the 
Black-Scholes equation and so can be directly applied to RSMR models. In fact Black- 
Scholes option pricing is used as a building block for more complex applications, for 



example Lieson Lei99l | values barrier options (options where the possibility to exercise 
depends on the underlying crossing a barrier) with jump risk using a Black-Scholes 
framework. 

Finally, the addition of regime switching a, also improves Fouque's option pricing 
method: 

• Firstly, Fouque calculates a from past historic data assuming no regime switching 
occurs (that is there is only one a value). This is because Fouque assumes 
stochastic volatility parameters do not change with time, which is not realistic 
over the long term. In fact it is worth noting that a current area of research 
or Fouq ue option pricing is addressing nonstationary stochastic volatility (see 



FPSOO 



• Secondly, it can be empirically shown that the expOU parameter values change 
with each regime; this is effectively the same as a changing with each regime 
i. This suggests we can improve Fouque's option pricing over the long term by 
introducing regime switching to model a. 

• Thirdly, a is fundamental to pricing C; in fact the most significant term in the 
expansion of C, Co, is the Black-Scholes option price with volatility a = a, 
therefore C is sensitive to changes in a. Specifically, the regime switching captures 
how d"j changes as the economy cycles through various economic phases and 
therefore should provide more accurate option prices. 

4. Numerical Experiment: Intra-Regime Option Pricing 
4.I. Calibration Procedure 

In this section we perform numerical experiments on pricing options on the S&P 
500 index. The purpose of the experiments is to firstly demonstrate we can obtain 
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intra-regime option prices for a general RSMR model, where its regime distributions 
are specified but its intra-regime dynamics are not. Secondly, we aim to compare the 
option pricing performance against Black-Scholes pricing and standard Fouque pricing. 
We demonstrate that by applying a regime specific CTj improves Fouque's option pricing, 
as opposed to applying one a which is currently proposed by Fouque. 

The experimental procedure was executed as follows. We fitted a linear plot to 
equation (11121) to extract a and h from the set of option prices for the chosen quote 
date; for "standard" Fouque option pricing we used only one a value for any quote 
period. For the Black-Scholes option pricing we set volatility a = a as we want 
to compare the option pricing performance of the Black-Scholes model for the same 
volatihty values. 

For our "regime based" Fouque option pricing, to calculate a and h in equation 
(11121) we applied the appropriate ai according to the regime i of the quote date. We 
identified the option quote date's regime using out of sample results from previous 
experimental results. To calculate ai it was shown in section 13.41 that is related to 
the return distribution's standard deviation for a given regime; these were calculated 
to be in regime 1 as 11.6%/year, regime 2 as 13.3%/ year. From previous experimental 
work the non-regime switching a was calculated to be a = 12.17%/ year. 

To compare the option pricing performance under Black-Scholes pricing, standard 
Fouque pricing and regime switching based Fouque pricing, we calculated the average 
percentage error of their prices against the empirically observed prices of the S&P 500 
options: 



where N is total number of option prices. Note that since regime switching has been 
proven to effectively model constant volatility over long periods, we expect Fouque's 
method to provide better results when using ai than just a. 

Numerical experiments performed with options data must be conducted differently 
compared to those on equities data. For a stock we can assign one price to it at a given 
point in time by taking the mid-point price of the bid and ask prices (as is generally 
practised in industry). A mid-point price is an acceptable indication of the stock's price 
since the bid-ask spread for stocks in general does not significantly vary over time. 

For an option assigning a single price for each point in time is more problematic than 
it is for stocks. Firstly, option prices are a function of a higher number of independent 
variables than stocks, so for a given point in time its quoted price will depend on the 
strike K and expiry T. Secondly, options are significantly affected by illiquidity effects 




I (observed option price) - (model option price) | 



observed option price 



(135) 



i=l 
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due to low trading volume. This leads to wide variation s in bi d-ask spreads, which in 
':u rn ar e a function of K and T (see for instance Pinder PinOSj, George and Longstaff 



GL93|). Consequently, there is no accepted method of assigning a single price to an 



option from bid and ask prices . 

Following Fouque |FPS00al | to avoid illiquidity effects we used highly traded options 
(S&P 500 options), used near at the money options with \K/ X — 1| < 3% and assigned 
the mid-point price of bid and ask prices for a given option price. We selected options 
data with expiries of at least three weeks since options close to expiry tend to be 
affected by more speculative and irrational effects. Under such effects Fouque states 
that it is doubtful if any diffusion model is useful for short expiries. Additionally, in 
section 13.41 it was proved that a becomes an increasingly better approximation of the 
stochastic volatility's constant volatility equivalent as time increases. Consequently, 
over shorter expiries a is not as good an approximation as the volatility process has 
not had sufficient time to frequently mean revert. 

4-2. State 1 Option Pricing Results 

We present the results of our option pricing numerical experiments in state 1. 



Table 1: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.14 


890 


43.3 


35.6 


45.5 


47.0 


0.14 


895 


39.9 


31.9 


40.7 


41.7 


0.14 


900 


36.6 


28.5 


36.4 


39.5 


0.14 


910 


30.4 


22.3 


28.8 


28.9 


0.14 


925 


22.4 


14.7 


20.0 


20.0 


0.14 


940 


15.7 


9.1 


13.7 


14.2 


0.14 


945 


13.8 


7.6 


12.2 


12.9 


0.21 


900 


43.2 


33.0 


42.0 


42.5 


0.21 


925 


29.1 


19.4 


26.3 


26.3 


0.38 


900 


54.4 


41.1 


52.1 


52.6 


0.38 


925 


40.8 


27.7 


37.0 


37.0 


0.63 


900 


66.3 


50.5 


64.2 


64.7 


0.63 


925 


52.9 


37.3 


49.5 


49.5 






Average Percentage 












Error 


28.8% 


6.9% 


6.4% 



Quote date: 29/4/03 
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Table 2: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.11 


975 


41.1 


36.0 


40.3 


41.7 


0.11 


980 


37.9 


32.2 


35.9 


36.7 


0.11 


995 


28.2 


22.3 


24.4 


24.6 


0.11 


1005 


22.9 


16.9 


18.3 


18.3 


0.11 


1025 


13.4 


8.8 


9.0 


9.3 


0.20 


975 


49.9 


41.6 


45.6 


46.5 


0.20 


980 


46.6 


38.2 


41.7 


42.4 


0.20 


985 


43.4 


35.0 


38.0 


38.5 


0.20 


995 


37.4 


28.9 


31.4 


31.5 


0.20 


1005 


31.9 


23.6 


25.5 


25.5 


0.20 


1025 


22.0 


15.0 


16.1 


16.3 


0.45 


1030 


19.9 


13.3 


14.1 


14.5 


0.45 


975 


64.9 


53.4 


57.6 


58.3 


0.45 


995 


52.9 


41.5 


44.8 


44.9 


0.45 


1005 


47.9 


36.3 


39.2 


39.2 


0.45 


1025 


38.0 


27.2 


29.5 


29.6 






Average Percentage 












Error 


23.0% 


16.5% 


15.6% 



Quote date: 7/7/03 
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Table 3: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.08 


995 


41.5 


32.5 


43.4 


45.1 


0.08 


1005 


34.6 


25.2 


33.1 


33.7 


0.08 


1010 


31.3 


22.0 


28.4 


29.0 


0.08 


1015 


28.3 


19.0 


24.8 


24.9 


0.08 


1025 


22.7 


13.7 


18.3 


18.3 


0.08 


1035 


17.8 


9.6 


13.1 


13.3 


0.08 


1050 


11.6 


5.1 


7.3 


8.3 


0.17 


995 


50.4 


38.8 


49.8 


50.9 


0.17 


1005 


44.0 


32.2 


41.4 


41.9 


0.17 


1025 


32.4 


21.2 


28.2 


28.2 


0.17 


1050 


20.4 


11.4 


17.0 


17.6 


0.25 


995 


56.4 


43.2 


55.0 


55.9 


0.25 


1005 


49.8 


36.9 


47.2 


47.6 


0.25 


1025 


38.3 


26.0 


34.5 


34.5 


0.25 


1050 


25.8 


15.7 


23.1 


23.5 






Average Percentage 












Error 


33.5% 


11.8% 


10.7% 



Quote date: 2/9/03 
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Table 4: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.12 


1005 


46.4 


34.2 


52.5 


53.9 


0.12 


1010 


42.9 


30.8 


46.8 


47.6 


0.12 


1015 


39.5 


27.5 


41.7 


42.2 


0.12 


1020 


36.5 


24.5 


37.1 


37.3 


0.12 


1025 


33.2 


21.7 


32.9 


33.0 


0.12 


1050 


19.5 


10.7 


17.0 


17.4 


0.12 


1060 


15.6 


7.8 


12.2 


13.1 


0.20 


1005 


53.1 


39.2 


58.3 


59.3 


0.20 


1015 


46.7 


32.8 


48.9 


49.3 


0.20 


1025 


40.5 


27.1 


40.9 


41.0 


0.20 


1035 


34.7 


22.1 


34.1 


34.1 


0.20 


1050 


27.4 


15.9 


25.7 


26.1 


0.20 


1055 


25.0 


14.1 


23.3 


23.8 


0.20 


1060 


22.6 


12.5 


21.1 


21.8 






Average Percentage 












Error 


36.0% 


7.2% 


6.9% 



Quote date: 6/10/03 
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The numerical experiments give option pricing results for state one, from four dif- 
ferent quote periods over a range of strikes and expiries. We see that regime based 
option pricing provides lower average percentage error compared to Black-Scholes or 
standard Fouque option pricing. In fact the average Black-Scholes option pricing error 
is 30.3% whereas Fouque's standard and regime based average option pricing errors are 
10.6% and 9.9% respectively. 

The average Black-Scholes option pricing error of 30.3% over different expiries and 
strikes validates Fouque's perturbation solution to stochastic volatility option pricing. 
Since the option price C is solved by the perturbation expansion equation (|5^ therefore 
C should be approximately equal to the leading order term Cq, which is the Black- 
Scholes option price. Additionally, it supports our claim that Cq can be used as an 
approximation for intra-regime option prices. 

Figure [H plots table[T]for options with expiry of 0.14 years, with the empirical option 
prices labelled "empirical". It is reassuring that all three pricing methods follow the 
same trend as empirical option prices, providing evidence that all three are viable 
option pricing methods. We observe that the Black-Scholes option prices consistently 
underprice over all strikes compared to the empirical prices, which has been frequently 
observed by many researchers. For instance, the Black-Scholes underpricing is most 
easily noticeable when empirical option prices exhibit volatility smiles. 

The Black-Scholes underpricing has been explained by researchers recognising that 
Black-Scholes option pricing assumes constant volatility and therefore no risk premium 
can be associated with any "volatility risk" . In stochastic volatility our risk premium 
associated with volatility is accounted for by the market price of volatility risk term, 
therefore it is possible to price options with volatility risk. In Fouque's perturbation 
solution to option pricing, Ci contains the market price of volatility risk, therefore 
suffers less from underpricing compared to the Black-Scholes equation. 
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4-3. State 2 Option Pricing Results 

We present the results of our option pricing numerical experiments in state two. 



Table 5: Option Pricing Results for SfcP 500 Call Option 



Rxnirv 


Strike 


Rmni rical 


Option Pricing Method 






Price 


rsl fi - Sir*!) nl pc; 


Perturbation Method 


1 Years 1 


fCcnts) 


fCents) 




St a n rl ^^rrl 


T~? po"irnp t;pn 


0.12 


1250 


75 fi 


4Q 7 


94 3 


92 8 


0.12 


1275 


59.4 


32 5 


58 4 


58 9 


0.12 


1300 


44.6 


19 4 


■^5 7 


■^5 6 


0.21 


1250 




Oc/ .O 


1 uu . o 




0.21 




73 Q 


49 7 


7^ 1 

I o. ± 


74 H 


0.21 


1285 


68 


OXJ.tJ 


65.4 


65.3 


0.21 


1300 


59.1 


29.2 


53 1 


53 1 


0.45 


1250 


122.6 


84 1 


143.7 


142.0 


0.45 


1275 


106.9 


67.6 


115.3 


114.7 


0.45 


1300 


91.8 


53.2 


93.1 


93.0 


0.72 


1250 


150.5 


107.0 


180.2 


178.0 


0.72 


1275 


135.0 


90.5 


152.1 


151.0 


0.72 


1300 


120.0 


75.6 


129.1 


128.7 


0.97 


1300 


145.4 


95.1 


160.0 


159.2 






Average Percentage 


39.4% 


11.2% 


10.6% 






Error 









Quote date: 2/1/01 
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Table 6: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.14 


1200 


66.9 


50.7 


82.5 


81.3 


0.14 


1210 


60.1 


43.5 


68.5 


67.8 


0.14 


1225 


50.6 


33.8 


51.7 


51.5 


0.14 


1250 


36.9 


20.8 


31.3 


31.3 


0.23 


1200 


80.2 


60.7 


93.4 


92.3 


0.23 


1225 


64.3 


44.3 


66.3 


66.1 


0.23 


1250 


50.3 


30.9 


46.7 


46.7 


0.40 


1200 


98.5 


75.6 


112.8 


111.7 


0.40 


1225 


83.0 


59.6 


87.9 


87.5 


0.40 


1250 


68.9 


45.8 


68.5 


68.4 


0.65 


1225 


105.4 


78.9 


112.8 


114.5 






Average Percentage 


32.1% 


9.6% 


9.0% 






Error 









Quote date: 26/4/01 
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Table 7: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.08 


1175 


50.9 


41.1 


62.3 


61.2 


0.08 


1200 


32.9 


23.6 


34.2 


34.1 


0.08 


1225 


20.5 


11.6 


17.5 


17.5 


0.17 


1175 


64.4 


50.5 


71.6 


70.8 


0.17 


1190 


55.0 


40.2 


56.7 


56.4 


0.17 


1200 


47.1 


34.1 


48.5 


46.3 


0.17 


1210 


41.3 


28.6 


41.3 


41.3 


0.17 


1225 


33.8 


21.5 


32.4 


32.3 


0.17 


1240 


26.9 


15.7 


25.3 


25.1 


0.42 


1175 


90.3 


70.1 


96.7 


96.0 






Average Percentage 


30.0% 


7.3% 


6.9% 






Error 









Quote date: 18/7/01 
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Table 8: Option Pricing Results for SfcP 500 Call Option 



Expiry 


Strike 


Empirical 


Option Pricing Method 






Price 


Black-Scholes 


Perturbation Method 


(Years) 


(Cents) 


(Cents) 


Pricing 


Standard 


Regime Based 


0.11 


1060 


63.4 


40.6 


80.6 


79.2 


0.11 


1070 


56.7 


33.2 


64.3 


63.7 


0.11 


1075 


53.3 


29.8 


57.7 


57.3 


0.11 


1080 


50.2 


26.6 


51.8 


51.6 


0.11 


1090 


44.1 


20.9 


42.1 


42.0 


0.11 


1100 


37.8 


15.9 


34.4 


34.4 


0.11 


1125 


26.3 


7.3 


21.7 


21.2 


0.18 


1100 


48.7 


22.3 


47.2 


47.2 


0.18 


1125 


36.0 


12.6 


34.6 


34.2 


0.28 


1160 


83.4 


52.3 


97.2 


96.4 


0.28 


1100 


59.0 


29.0 


60.5 


60.5 


0.28 


1125 


46.1 


18.6 


47.5 


47.3 






Average Percentage 












Error 


51.5% 


9.4% 


9.1% 



Quote date: 10/9/01 
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Figure 2: Graph of Table [5] for Optfons with Expiry of 0.11 Years 
Call Option Prices for S&P 500 (Expiry 0.1 1 Years) 
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The numerical experiments give option pricing results for state two, from four 
different quote periods over a range of strikes and expiries. As in state one regime 
based option pricing provides lower average percentage error compared to Black-Scholes 
or standard Fouque option pricing. In fact the average Black-Scholes option pricing 
error is 38.2% whereas Fouque's standard and regime based average option pricing 
errors are 9.4% and 8.9% respectively. Therefore the pricing error under Black-Scholes 
option pricing is higher in state two than in state one, which is consistent with empirical 
observations as volatility smiles increase during "down states" (state two) and therefore 
option prices deviate away from a constant volatility assumption. 

Figure [2] plots table for options with expiry of 0.11 years with the empirical 
option prices labelled "empirical". As in state one, all three pricing methods follow 
the same trend as empirical option prices, providing evidence that all three are viable 
option pricing methods in state two. We also note the Black-scholes option prices are 
consistently underpricing over all strikes, as in state one, for the same reasons as in 
state one -Black-Scholes option pricing does not take into account "volatility risk" . 

From figure [2] we observe that both Fouque based option pricing methods exhibit a 
more pronounced volatility smiling effect compared to figure [H Since we price options 
upto the first correction \feC\ and Cq does not take into account any volatility smiling 
effects, we can attribute the smiling effect to \feC\. The term y^Ci can increase 
volatility smiling through an increase in the market price of volatility risk, which is 
likely because state two models the "down" state and we expect risk associated with 
volatility to be higher. Additionally, we notice that \fiC\ contains p, the correlation 
between the volatility's Wiener process W2 and which may increase in "down" 
states compared to "up" states. This is consistent with empirical observations since 
volatility tends to increase during a down state. 

5. Conclusions 

It can be seen from our literature review of volatility models that the development 
of volatility models has progressed in a logical order to address key shortcomings of pre- 
vious models. Time dependent models addressed option prices varying with expiration 
dates, local volatility also addressed volatility smiles and the leverage effect, whereas 
stochastic volatility could incorporate all the effects captured by local volatility and a 
range of other empirical effects e.g. greater variability in observed volatility. However 
the trade-off associated with improved volatility modelling has been loss of analytical 
tractability. 

In conclusion we have shown that Fouque's option pricing method provides signif- 
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icantly better option pricing accuracy compared to Black-Scholes option pricing over 
a variety of strikes and expiries. We have provided option prices and have shown that 
taking into account regime specific a^, rather than using one a, improves option pricing 
accuracy. 
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